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FREQUENCY SHIFTS INDUCED IN LASER PULSES 
BY PLASMA DENSITY VARIATIONS 


The propagation of electromagnetic radiation through plasmas is a problem of gen¬ 
eral interest with a wide variety of applications ranging from communications 1 to laser 
driven particle accelerators . 2 ’ 3 For example, the results of plasma simulation studies have 
recently been reported 4 ’ 5 which suggest two possible methods by which the frequency 
of an electromagnetic wave may be upsliifted. In the first method, the plasma density 
through which the EM wave is propagating is suddenly increased in time , 4 while the sec¬ 
ond method utilizes the interaction of a plasma wave (having a phase velocity near the 
speed of light ) 2 ’ 3 ’ 6 ’ 7 with a short EM pulse . 5 Phenomena such as these, which result from 
variations in the plasma density, may offer a way of tuning the radiation from a laser or, 
alternatively, they may describe the distortion of radio signals in the ionosphere. 

In the following, a general analytic theory is developed from first principles which 
describes how a EM radiation pulse is affected by variations in the plasma density. Specif¬ 
ically, the ID wave equation is used to derive expressions for the shifts induced in the 
frequency and ia the wavenumber by arbitrary plasma density variations (both in space 
and in time). In particular, for variations in the plasma density as a lunction of time, it is 
shown that the wavenumber of the EM wave remains constant while shifts are induced in 
the frequency. (This is to be contrasted to the case of spatial plasma density variations, 
for which the frequency remains constant and the wavenumber shifts.) For the case of a 
short radiation pulse interacting with a plasma wave with finite phase velocity, it is shown 
that shifts are induced in both the wavenumber and in the frequency. Maximum frequency 
shifts may be obtained when the phase velocity of the plasma wave is equal to the group 
velocity of the radiation pulse evaluated at the ambient plasma density. Furthermore, it 
is shown that for small frequency shifts, the amplitude of the vector potential maintains 
its initial profile while propagating forward at the ambient group velocity. 

The ID wave equation for the normalized vector potential a = eA/(mc 2 ) of the radi¬ 
ation field is given by 

a(z,<) = (i) 

where A ; 2 = Ar 2 „n(z, t)/(‘y(z, t)n 0 ). Here 7 (z,<) is the relativistic factor associated with the 
motion of the plasma electrons, n(z,<) is the plasma electron density and A ; 2 0 = u> 2 0 /c 2 , 
where u > p0 is the electron plasma frequency in the ambient density n 0 . In deriving the 
above equation, conservation of canonical momentum was used, p± — cA±/c, which gives 
a trans.-erse plasma current (in the fluid limit) of J± = —enp±/(mj). Throughout the 
Manuscript approved September 8, 1989. 
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following it will be assumed that a 2 << 1 such that k 2 (z,t) is independent of a(z,t), i.e., 
the effects of the radiation field on the plasma wave will be ignored. 

To solve the above wave equation, it is helpful to write a(z,t) = 6(2,t)exp(ifco~ — 
iu>o<), where b(z,i) is the radiation envelope and where u>o and k 0 are the frequency and 
wavenumber of the radiation in the ambient plasma (in the absence of a plasma wave) 
which satisfy the dispersion relation u)$ = c 2 k% + c 2 k 2 0 . Furthermore, it is convenient to 
introduce a change of variables ( = z — v g t and r = t, where v g = c 2 k 0 /u 0 is the group 
velocity of the radiation in the ambient plasma. The wave equation is then given by 


( 


.td 0 0 
~!h + 


v g d 2 1 d 2 

^dCdT + ^ g dl 2 


i d 2 \ 

c 2 dr 2 ) 


6 = 




( 2 ) 


•.vhere bk 2 = k\- * 2 0 and I/7 2 = i - v]/c 2 . 

It will now be assumed that the radiation envelope 6(C,r) is slowly varying compared 
to the radiation frequency u> 0 , that is \db/dr\ << |u> 0 6| and \db/d(,\ << |u> 0 6/c|. Assuming 
this, the second order derivatives in Eq. ( 2 ) may be neglected. Equation ( 2 ) may then be 
solved giving 

6«,t) = &o«)exp . ( 3 ) 

The condition \db/dr\ << (u^o^l implies \c 2 bk 2 /( 2 u)l )| << 1. For a plasma density pertur¬ 
bation 6 n, this condition implies |u7 2 0 6n/(2u^n 0 )| << 1 . The condition \db/d(,\ « |u? 0 6/c| 
implies | / Q r dT'{dbk 2 /d£)c 3 /( 2 ul)\ << 1 as well as |(d6 0 /dC) c /(&o u, o)| << 1 . The first of 
these inequalities generally implies that the wavenumber shift must be small compared to 
<*>o/c, while the second inequality indicates that the initial envelope 60(C) must be slowly 
varying compared to u>o/c. 


The total phase $(C,t) of the radiation field may be identified by writing a(C,T) = 
|a((,T)| expi$(C,T). It is then possible to examine the evolution of the frequency as 
well as the wavenumber of the radiation through the definitions u{(,r) = -d$/dt = 
—(d/dr — Vgd/dQfy and A:(£,r) = d$/dz = d$/d(. Using the above solution for 6 (C,t) 
gives 

U,(< ’ T) = ““ + ^^ (< ’ T) - ^ l dr T ( sk l«' T,) ' < 4 > 


KCr)=k,, ~h,L dT 'r( Sk ^ y) - (5) 

The above equations are valid for arbitrary variations 6fc 2 (C,r). (Recall, k 2 = k 2 0 n/(n 07) 
and the effects of 7 become important for nonlinear relativistic plasma waves. 7 ) Provided 
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ft(£,r) remains slowly varying compared to wo, the amplitude of the normalized vector 
potential |a| does not change. That is, Eq. (3) indicates that |a| = IMOI and, hence, the 
initial envelope of the vector potential is simply convected forward at the group velocity 


To illustrate the above theory, consider a plasma density variation which is a function 
only of space. For example, consider a radiation pulse entering a plasma (w 2 p /*l « 1 ) 
from vacuum with a plasma density profile 8k 2 = 8k 2 (z) for z > 0 and equal to zero 
for z < 0. Assume that at t = 0, the radiation pulse extends from —L < z < 0, where 
a > 0 L/c » 1. Equations (3)-(5) indicate, that as the pulse propagates, the frequency and 
wavenumber evolve according to *(z,t) = wo and k(z,t) = fc 0 — c 2 8k p (z)/(2v g u) 0 )- This 
is in agreement with the well-known result 1 that as radiation propagates into a plasma 
with spatial density variations, the frequency remains constant whereas the wavenumber 
changes such that the dispersion relation w 2 = c 2 (k 2 + k 2 ) remains satisfied. 

On the other hand, consider a plasma density variation which is a function only of time. 
For example, consider a radiation pulse (of length L >> c/w 0 ) propagating through a long, 
uniform plasma column (where * 2 (*l << 1 ) in which the density is temporally changing, 
6k 2 = 8k 2 (t). Equations (3)-(5) indicate, that as the pulse propagates, the frequency and 
wavenumber evolve according to u>(z,t) = w 0 + c 2 8k 2 (z)/( 2w 0 ) and fc(z,<) = k 0 . This is 
in agreement with the simulations of Ref. 4 which indicate that as radiation propagates 
through a plasma with temporal density variations, the wavenumber remains constant 
whereas the frequency changes such that the dispersion relation u> 2 = c 2 (k 2 + k 2 ) remains 
satisfied. 

Consider a plasma density variation which is a function both of time and of space. 
For example, consider the case in which a plasma wave (with phase velocity near c ) 2 ' 3,6,7 
is used to upshift the frequency of a laser pulse, as suggested by the simulations of Ref. 5. 
Assuming a plasma wave with a phase velocity v p such that the plasma wave quantities 
are functions of only z - v p i implies 8k 2 (C,r) = 6k 2 (( - Avr), where 8k 2 (( - A vt) has 
the form of a periodic oscillation and Aw = v p - v g . Defining the shift in frequency and in 
wavenumber as Aw = w(£,r) — w((,0) and A k = fc(£,r) — fc(CiO) gives 

Au,= 2 Sj 1< *' (< “ A '’ T) "^ (0! ’ (6) 

At = " a ” t) " ’ (T » 

The above equations indicate that the frequency shift will be maximum for the case 
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v p = Vg. For this case Aw = — (v g rc 2 /2u> 0 )d6k p /d( and A k — —(c 2 r/2uo)d6k p /d(. Hence, 
Aw is a linear function of v g r, the distance that the pulse has propagated through the 
plasma. In such a way |Aw| will increase until it becomes sufficiently large so that the 
approximation that b((,r) is slowly varying is no longer valid. For 6(C,r) to be slowly 
varying requires |cAw/(v ff w 0 )| 2 << 1. Consider a plasma wave with v p = v g such that 
6k 2 (() = k p0 6n(()/no, where 6n(() is the plasma wave density perturbation For 6n(() = 


Sno sinfcpoC? then 


Aw 


o Sn 0 v g ' 


- = — ir —\ -cos k p0 (. 

Wo Wq Tlo A p o 


( 8 ) 


Hence, a positive frequency shift Aw > 0 requires the laser pulse to be positioned prop¬ 
erly in the phase of the plasma wave such that dSn/dC, < 0. This is illustrated in Fig. 
1. A convenient method for producing such a plasma wave may be the laser Wakefield 
accelerator, 3 in which a driving laser pulse of frequency w 0 is used to generate a plasma 
wave with v v = v a . 


For the case v p ^ v g , |Aw| is no longer a linearly increasing function of cr. In fact, 
Eq. (6) indicates that Aw will oscillate as a function of cr. Consider a plasma wave of 
the form Sk 2 ((,r) = k' 0 \Sn/n 0 \cos k p0 (( - Avr). Assuming the laser pulse (with a pulse 
length L « A p0 = 2ir/k p0 ) is initially centered about Co = 0, Eq. (6) indicates that 
the maximum frequency shift is given by |Aw m | = \u> p0 v p 6no /(u>on 0 Av)\ and occurs when 
cr = |Ap 0 /(2Av)|. Furthermore, notice that Aw m may either be positive or negative, 
depending on the sign of Av (for the present example, Aw m > 0 for v p < v g ). Again, the 
assumption that b((,r) is slowly varying implies |cAw m /(2vpW 0 )| 2 << 1. 


Physically, the frequency shifts induced in a radiation pulse by a plasma wave may 
be understood as follows. A plasma wave gives rise to variations in the plasma parameter 
k p {C,,r). These variations lead to variations in the “local” phase velocity of the laser 
pulse. Heuristically, the “local” dispersion relation for the radiation field is given by 
w 2 «,r) = c 2 k 2 ((,r) + c 2 k 2 ((,r). For small cr, this gives v 2 {C)/c 2 = 1 + kl o n(()/{k^n 0 ). 
For example, the local phase velocity near the front of the laser pulse (( = Co) will be less 
than the local phase velocity near the back of the pulse (C = Cl) provided n(Co) < M(l)- 
Hence, the individual phase peaks in the pulse a(z,t) = |a|exp(iA;z - iwf) may move 
relative to one another (i.e., closer together for the present example). In such a way both 
the frequency and the wavenumber of the radiation pulse will change, as is given by Eqs. 
(4) and (5). 

It should be pointed out that due to the local nature of the frequency shift (the 
dependence of Aw on C), a laser pulse with a finite pulse length L will develop a spread 
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in frequency shifts. That is, the frequency shift at the front of the pulse may be different 
than the frequency shift at the back of the pulse. For example, consider a plasma wave 
with t>p = v g of the form 8k 2 = fc p0 (6no/n 0 )sin fc p0 C and a laser pulse centered about £ = n 
with L < Apo/2 (see Fig. 1). The difference in the frequency shift at the center of the pulse 
with a point AC away is given by |(Au»(C) - Au>(( + AC))/Au;(C)| = 1 - cos AC- This 
spread in frequency shifts may be quite significant. 

Equations (4) and (5) provide convenient analytic expressions for the frequency and 
wavenumber shifts resulting from variations in the plasma parameter k 2 . However, these 
expressions are only valid for a slowly varying envelope 6(£,r), which implies that the 
frequency shifts must be small, |Aa>/u;o| 2 << 1. To calculate the behavior of the vector 
potential for large shifts, |Au>/u>o| > 1, it is necessary to retain the d 2 /(d(dr) derivative 
in the wave equation for b( C,t). Keeping this derivative, the wave equation can formally 
be solved for the case 6k 2 (C,r) = 6k 2 (£), giving 

b (C» T ) = ^ ^C'-FXOexpftuVC' - C)/^]^o[o:(C\C, T )]> (9) 

a«’,(,T)= ^2cV dCSkKO/v, 

where J 0 is the Bessel function of order zero and F(() = (d/d( + iu> 0 /v s )i> 0 ((). In solving 
for &(C,r), it has been assumed f»(C = 0) = 0, i.e., C = 0 is chosen to be in front of 
the pulse. Analysis of the above equation indicates that now both the amplitude |fr(C> r )| 
and the phase evolve as a function of r. Determination of the precise behavior of b((,r) 
requires further computation. 



In the analysis of Ref. 5, heuristic arguments based on Lorentz transformations of the 
dispersion relation are used to determine an expression for dw/dz. The above analytic 
theory may be used to derive an equation for dw/dz by operating on Eq. (4) with djdz = 
d/d(. For the case of a plasma wave with phase velocity v p , Eq. (4) gives 


du> 

dz 


2u> 0 Av 


- v P t) - v g—6k‘{z - 


( 10 ) 


which for v p = v g reduces to the expression given in Ref. 5 when 6k 2 = k 2 0 6n/n 0 . Hence, 
the expression given in Ref. 5 has been shown to be valid provided i) v p = v g , ii) v p is 
independent of z and t, iii) relativistic effects associated with 7 (z,t) are neglected and iv) 
the frequency shift remains small, |Au>/u > 0 | 2 << 1 . Whether or not the arguments of Ref. 
5 are valid for large shifts |Au>/u>o| > 1 remains uncertain. 
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The analytic theory presented above describes how a laser pulse becomes modified due 
to variations in the background plasma. Specifically, arbitrary variations in the plasma 
parameter kp((,r) leads to shifts in both the frequency and wavenumbers indicated by 
Eqs. (4) and (5). Provided these shifts are sufficiently small, the amplitude of the vector 
potential remains unchanged and propagates forward at the group velocity v g . The above 
theory indicates that changes in the plasma density as a function of time lead to shifts in the 
radiation frequency. In particular, the possibility of upshifting the frequency of the laser 
pulse using a plasma wave has been examined and it is found that maximum frequency 
shifts result when v p = v g . For this case, Aw increases linearly with the propagation 
distance v g r and positive frequency shifts require phasing the laser pulse such that it 
is centered at a position of decreasing density or, more precisely, at a position where 
d6kp/d( < 0. The above theory is valid, in general, for small frequency shifts, |Aw/w 0 | 2 < < 
1. However, Eq. (8) implies that large shifts of order |Aw| ~ wo require propagation 
distances of order cr ~ |<Jno/no| -1 w^/wp 0 . Since, typically (for a KrF laser), w^/w 2 0 >> 1, 
this propagation distance can be quite large. The detailed examination of large shifts 
| Aw | ~ w 0 requires further evaluation of Eq. (9). 

This work was supported by the Department of Energy and the Office of Naval Re¬ 
search. 
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